Finite- Temperature Phase Diagram of Hard-Core Bosons in Two Dimensions 
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We determine the finite-temperature phase diagram of the square lattice hard-core boson Hubbard 
model with nearest neighbor repulsion using quantum Monte Carlo simulations. This model is 
equivalent to an anisotropic spin-1/2 XX Z model in a magnetic field. We present the rich phase 
diagram with a first order transition between solid and superfiuid phase, instead of a previously 
conjectured supersolid and a tricritical end point to phase separation. Unusual reentrant behavior 
with ordering upon increasing the temperature is found, similar to the Pomeranchuk effect in '^He. 
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A nearly universal feature of strongly correlated sys- 
tems is a phase transition between a correlation-induced 
insulating phase with localized charge carriers, to an 
itinerant phase. High temperature superconductors [Q, 
manganites [|| and the controversial two dimensional 
(2D) "metal-insulator-transition" ^ are just a few ex- 
amples of this phenomenon. The 2D hard-core boson 
Hubbard model provides the simplest example of such a 
transition from a correlation-induced charge density wave 
insulator near half filling to a superfiuid (SF). It is a pro- 
totypical model for preformed Cooper pairs of spin 
flops in anisotropic quantum magnets ||^,^ , of SF Helium 
films and of supersolids |^,^ . 

In simulations on this model, which does not suffer 
from the negative sign problem of fermionic simulations, 
we can investigate some of the pertinent questions about 
such phase transitions: what is the order of the quan- 
tum phase transitions in the ground state and the finite 
temperature phase transitions? Are there special points 
with dynamically enhanced symmetry [lO|]? Can there be 
coexistence of two types of order (such as a supersolid - 
coexisting solid and superfiuid order)? Answers to these 
questions provide insight also for the other problems al- 
luded to above. 

The Hamiltonian of the hard-core boson Hubbard 
model we study is 



(iJ> 
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where O; (ai) is the creation (annihilation) operator for 
hard-core bosons, ni ~ OjAi the number operator, V the 
nearest neighbor Coulomb repulsion and fi the chemi- 
cal potential. This model is equivalent to an anisotropic 
spin-1/2 XX Z model with = V and \Jxy \ = 2i in a 
magnetic field h — 2V — /i. The zero field (and zero mag- 
netization = 0) case of the spin model corresponds to 
the half filled bosonic model (density p = {rriz) = 1/2) at 
ji = 2V . Throughout this Letter we will use the bosonic 
language, and refer to the corresponding quantities in the 
spin model where appropriate. Due to the absence of ef- 
ficient Monte Carlo algorithms for classical magnets in a 
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FIG. 1. Ground-state phase diagram of the hard-core bo- 
son Hubbard model. The dashed line indicates the cut along 
which we calculated the finite-temperature phase diagram 
shown in Fig. ^. 



magnetic field there are still many open questions even in 
the classical version of this model, which was only studied 
by a local update method ||ll[]. 

In Fig. H we show the ground-state phase diagram 
flplJT^. For dominating chemical potential /i the sys- 
tem is in a band insulating state (p = and p = 1 re- 
spectively), while it shows staggered checkerboard charge 
order (p = 1/2) for dominating repulsion V. These 
solid phases are separated from each other by a SF. Ear- 
lier indications for a region of supersolid phase between 
checkerboard solid and SF phase turned out to be due to 
phase separation at this transition which is of first order 
atT = |,Ul]]. 

All of these phases extend to finite temperatures. On 
the strong repulsion side the hard-core boson Hubbard 
model is equivalent to an antiferromagnetic Ising model 
at i = and the insulating behavior extends up to a 
finite-temperature phase transition of the Ising univer- 
sality class (Tc « 0.567F at half filling). In case of zero 
repulsion the system is equivalent to an XY model in 
a magnetic field and the SF extends up to a transition 
of Kosterlitz Thouless (KT) type at finite temperature. 
The point V = 2t and p = 2V is also of special interest. 
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it corresponds to the Heisenberg point in the spin model 
and has an enlarged 0(3) symmetry, with long range or- 
der only at zero temperature in 2D [T^ ]. 

We focus our attention on the finite-temperature phase 
diagram in the vicinity of the first order quantum phase 
transition separating the SF and solid ground states away 
from the Heisenberg point. It is instructive to first con- 
sider the three dimensional (3D) version of this model, 
with a similar ground-state phase diagram. There, Fisher 
and Nelson ||T^ determined that the first order line ter- 
minates at a finite-temperature bicritical point at which 
the transition temperatures of the solid and SF phases 
meet. This bicritical point has a dynamically enhanced 
0(3) symmetry - the symmetry of the Hamiltonian is 
only 0(2). 

This feature alone dictates that despite the simi- 
lar ground-state phase diagrams, the finite-temperature 
phase diagrams have to be very different in 2D, as an 
0(3) critical point in 2D can exist only at zero temper- 
ature . Even with the absence of a supersolid in the 
ground state several alternative scenarios are still plau- 
sible: (i) the bicritical point could be pulled down to 
zero temperature, and the transition temperatures Tkt 
and Tc approach zero temperature without a direct phase 
transition between SF and solid at finite temperatures, 
(ii) a supersolid phase could exist at finite temperatures, 
although it is absent in the ground state, or (iii) the phase 
transition lines meet at finite temperature, but without 
dynamically enhanced symmetries. 

In order to decide between these possibilities and deter- 
mine the finite-temperature phase diagram we use quan- 
tum Monte Carlo simulations, employing a variant [^jl5| 
of the worm algorithm in the stochastic series ex- 
pansion (SSE) representation The simulations were 
performed in the grand canonical ensemble using lattices 
up to size 96 x 96 and do not suffer from any systematic 
error apart from finite size effects. In contrast to the loop 
algorithm ||l7|] the worm algorithm, while not as efficient 
at half filling (a spin model in zero magnetic field), re- 
mains effective away from half filling (in a finite field), 
where the loop algorithm slows down exponentially. 

From now on wc will restrict our discussion to = it, 
which corresponds to the dashed line in Fig. |l|. We 
present the finite-temperature phase diagram for this rep- 
resentative value in Fig. |2[ Due to particle-hole symme- 
try it is sufficient to show densities up to 1 /2 and chemical 
potentials up to /i = 2V = 6i. We find that scenario (iii) 
is realized: the first order phase transition extends to a 
tricritical point (D) at around « 0.4i, where it meets 
with the second order melting transition of the solid. The 
KT transition temperature Tkt of the SF meets the first 
order line at a lower temperature (B) ~ in contrast to the 
3D case, where all the three lines meet at the same bi- 
critical point (B=D) with enhanced 0(3) symmetry. We 
want to emphasize one unusual feature before discussing 
in detail how the phase diagram was determined: there is 
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FIG. 2. Phase diagrams a.t V = 2)1 as function of a) /i and 
T and b) p and T. The arrow in the inset of the upper panel 
denotes an upper bound for the critical /j. sd, T — 0.5t. The 
error bars of the triangles are strict bounds for the KT tran- 
sition. The right and upper axis show the parameters in the 
spin language. The solid circle and square denote qualita- 
tively where we expect the critical points to be. Lines are 
guides to the eye. 



reentrant behavior in the vicinity of the tricritical point 
(D) and the system melts upon decreasing the tempera- 
ture at constant chemical potential. 

The first order phase transition is characterized by a 
finite jump Ap in the density. As discussed in Ref. [p^ , 
where we present the ground-state phase diagram, two 
peaks can be observed in the density distribution at the 
critical chemical potential. Since a two-peak structure 
can also be seen at a second order transition in a finite 
system, a careful finite size scaling analysis is needed. 
Figure |3^ shows the finite size scaling of the densities of 
the two coexisting phases. For T < 0.39t the jump Ap 
remains finite under extrapolation, thus confirming the 
first order nature up to this temperature (line from C to 
D in Fig. 0), while it scales to zero for T > OAbt, an in- 
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FIG. 3. a) Finite size scaling of the coexistence densities at 
different temperatures, demonstrating the first order nature 
of the phase transition at T < 0.39t. b) Finite size scaling 
of the SF density ps at different temperatures, showing that 
the first order transition is between a solid and SF at low 
temperatures, and between a solid and normal fluid at higher 
temperatures. The dashed line shows the value of ps/(Tm) 
at a KT transition. All simulations are done at V = St. 



dication for a second order transition (line from D to E). 
Additionally we have checked the finite size scaling of the 
fourth order cumulant ratio C4 = 1 — (mf )/3(m^)^, of the 
staggered charge order parameter in the solid phase. 
Consistency with the expected size independent univer- 
sal value of the Ising universality class C4 = 0.6106900(1) 
p8[ at the transition confirms the second order nature. 
We have thus determined that the first order transition 
extends up to a critical point (D) at a temperature larger 
than 0.39t and smaller than 0.45t. In the next paragraph 
we show that the critical point where the KT line meets 
the first order line (point B) is below 0.39t. Therefore the 
end point D of the first order line is a iricritical point. 

The KT transition temperatures Tkt (line from A to B 
in Fig. 1^) can be determined using the SF number density 
Ps = {mT/2){W2 + W^) m. Here and Wy are the 
winding numbers of the world lines in x and y direction 
respectively, and m = l/(2t) is the mass of a boson. For 
system size L — 00, ps has a universal value at Tkt, 
Ps{Tkt) — (2/7r)mTifr, while it is always larger than 
this universal value for T < Tkt- At T = Tkt, Ps jumps 
from the universal value to zero and we have Ps{T) = 
for all T > Tkt- However, for finite L, ps is nonzero at 
all T and the finite size corrections of ps at T = Tkt are 
given by |§ p,7r = 2mTKTil + {2\og[L/ Lo{Tkt)]}-^)- 
Although these finite size corrections are known to be 
notoriously difficult, we can obtain strict upper and lower 
bounds for Tkt by plotting l/^nps/mT — 2) — log(i) 
as function of logL. As L is increased, this quantity 
converges to the constant (— logLo) at Tkt and diverges 
to ±00 for T > Tkt and T < Tkt, respectively [pi. 
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FIG. 4. a) Finite size scaling of the critical chemical poten- 
tial fic st T — 0.39t. The horizontal line is an upper bound 
for the critical chemical potential at T = 0.5t. b) Finite size 
scaling of the fourth order cumulant C4 at T = 0.5t for differ- 
ent chemical potentials. The horizontal line is the universal 
value of C4 at the Ising critical point. All simulations are 
done at V = 3t. 



The error bars we obtain using this method are reliable 
and small enough for our purposes. 

To investigate how the KT transition line meets the 
phase separation line we perform canonical measure- 
ments of the SF density ps at the density p/ of the fluid 
on this coexistence line. In Fig. we show the finite 
size scaling of ps{pf)- At low temperatures T = 0.3t ps 
of the SF at the coexistence line is well above the univer- 
sal value for a KT transition. Therefore the transition 
between SF and solid at temperatures T < 0.3t does not 
include a KT transition. The system remains SF right 
up to a direct first order phase transition from the SF 
to the solid. At T = 0.39t on the other hand, ps as a 
function of the system size L drops below the universal 
value and should scale to zero according to |^l|]. There 
is thus a first order transition from solid to normal fluid 
at r = 0.39t and the KT line has to meet the first order 
line at a finite temperature between 0.3t and 0.39i. 

As can be seen in Fig. ||, the checkerboard solid is 
melting in a first order transition to a SF upon de- 
creasing the temperature around p ~ 2.35i. Even more 
special is the melting to a normal fluid in the vicinity 
of the tricritical point (D) both upon increasing or de- 
creasing the temperature. To proof this reentrant be- 
havior, we first determine the critical chemical potential 
Pc(T = 0.39i) of the first order transition at T = 0.39t. 
We then show that pc{T = 0.5t) < p^T = 0.39i) On 
a finite system we define pc as the chemical potential 
where there is an equal probability to be in the solid and 
fiuid phase (equal peak areas in the density probabili- 
ties). Performing a finite size extrapolation (Fig. ^) we 
obtain /L(c(0.39t) = 2.3015(10)i. Next we calculate the 
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upper bound for the critical chemical potential /ic(T) at 
T = QM. Observing that the fourth order cumulant 
ratio C4 scales above the universal critical value C4 at 
^l = 2.295t (Fig. |b), we thus have Atc(T = 0.5t) < 2.295t. 
As ^c(r — Q.bt) < fic{T = 0.39t) we have pinpointed 
reentrant behavior where the solid melts to a normal fluid 
both when increasing and decreasing the temperature. 

To summarize, we have investigated the finite- 
temperature phase diagram of the hard-core boson Hub- 
bard model with nearest neighbor repulsion using quan- 
tum Monte Carlo simulations. The first order phase tran- 
sition in the ground state between SF and solid phases 
extends to finite temperatures and turns into a second or- 
der melting line of the solid phase at a finite-temperature 
tricritical point. As the KT line meets the first order line 
at a lower temperature, we have a qualitatively different 
phase diagram than in 3D [0 , where there is a bicritical 
point with dynamically enhanced 0(3) symmetry. 

We expect this phase diagram, obtained for a represen- 
tative value otV — 3t, to be valid also for other couplings 
V > 2t. On approaching the Heisenberg point V ^ 2t the 
temperature scale of all the transitions drops to zero. In 
the other limit t/V — > 0, we approach the antiferromag- 
netic Ising model. The KT transition temperatures Tkt 
and the tricritical point should scale with t, since in 
the Ising model the phase transition is first order at zero 
temperature, but second order at any finite temperature 
lUl, consistent with T3/V - t/V 0. 

The melting line of the solid shows several unusual 
features. As a function of chemical potential there are 
three different phase transitions: at low temperatures the 
melting transition is first order into a SF. At intermediate 
temperatures the melting transition is still first order, but 
to a normal fluid, while it is a second order transition to 
a normal fluid at even higher temperatures. 

There is reentrant behavior in the vicinity of the tri- 
critical point: as a function of increasing temperature 
there are three phase transitions: first the SF turns into 
a normal fluid. Upon increasing the temperature further 
this normal fluid solidifies, and finally melts again. The 
reason for this reentrant behavior is conjectured to be 
a large entropy of the weakly interacting vacancies (or 
interstitials at p > 1/2) in the solid, larger than the en- 
tropy of the normal fluid. Such an ordering transition 
upon increasing temperature is rare but has been ob- 
served before, e.g. in the bcc Ising antiferromagnet p3| . 

Finally, we want to mention the striking qualitative 
similarity of the presented phase diagram to those of 
fermionic "^He and bosonic ^He in two dimensions. In 
■^He there is also a reentrant melting line (the Pomer- 
anchuk effect [Q) and a SF phase at low temperatures 
psf . In '^He however the large entropy of the solid is due 
to the nuclear spin degrees of freedom and the effect even 
more pronounced. In ^He films, the solid phase melts to 
a SF in a small range for the pressure upon decreasing the 
temperature pw, like in our phase diagram where there 



is a direct transition form SF to solid for constant values 
of fj, around /x ~ 2.35t. Note that in our phase diagram 
we have pinned down the even more subtle effect of a 
reentrant melting line between solid and normal fluid. 
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